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On a fractional linear birth-death process 

ENZO ORSINGHER* and FEDERICO POLITO** 

^Department of Statistics, Probability and Applied Statistics, Sapienza University of Rome, pi. A. 
Mora 5, 00185 Rome, Italy. E-mail: enzo.orsingher@uniromal.it; federico.polito@uniromal.it 

In this paper, we introduce and examine a fractional linear birth-death process N^it), t > 0, 
whose fractionahty is obtained by replacing the time derivative with a fractional derivative in the 
system of difference-differential equations governing the state probabilities Pfe(f), t>Q, k > 0. 
We present a subordination relationship connecting Nv{t), t>0, with the classical birth-death 
process N{t), t> 0, hy means of the time process T2v{t), t>0, whose distribution is related to 
a time-fractional diffusion equation. 

We obtain explicit formulas for the extinction probability Po{t) and the state probabilities 
p%{t),t> 0, k> 1, in the three relevant cases A>/i, A</i, A = /i (where A and fj, are, respectively, 
the birth and death rates) and discuss their behaviour in specific situations. We highlight the 
connection of the fractional linear birth-death process with the fractional pure birth process. 
Finally, the mean values ¥.N,^{t) and Var Af^(f) are derived and analyzed. 

Keywords: extinction probabilities; fractional derivatives; fractional diffusion equations; 
generalized birth-death process; iterated Brownian motion; Mittag-Leffler functions 



1. Introduction 

In a previous paper [8], we constructed a fractional version of the pure birth process 
Afi,{t),t > (both in the general and in the linear case denoted here as M^{t),t > 0), 
by considering the fractional equations governing their distributions. In this work, we 
examine the linear birth-death process N,j{t), t>0, where the state probabilities 

pUt)^Pr{N.{t)^mAO) = l} (1-1) 
are assumed to satisfy the fractional difFerence-difFerential equations 

= - (A + ^^)kpk{t) + X(k - l)pk-i{t) 

(1.2) 

+ ^i{k + l)pk+i{t), k>l,0<iy<l. 
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The fractional operator appearing in (1.2) is defined as 
dVW 1 (d/d,s)/(g) ^ 

— ; = — -, r / — -, ^ ds, [)<v<l, n o\ 

At'' T{l-v)Jo {t-sy ' (1-3) 

fit), v = l. 

The derivative (1.3) is usually ealled a Caputo or Dzherbashyan-Caputo fractional deriva- 
tive and differs from the classical Riemann-Liouville derivative by exchanging the inte- 
gral and derivative operators (see [9]). An advantage of Caputo over Riemann-Liouville is 
that Caputo does not require fractional-order derivatives in the initial conditions, which 
is good for practical purposes. The positive parameters A and are, respectively, the 
birth and death rates. 

The exact distribution of the linear birth-death process reads (see [1], page 91, [4], 
page 454) 

Plit) = (A - A^)^c-(^-^)^ ^^;^_-;;_^^_^^/^i^^ , > l,t > 0,^ ^ A. (1.4) 

When A = /i, the distribution (1.4) is much simpler and takes the form 

pl{t)= , , t>0,fc>l. (1.5) 

' (l + Ai)'=+i' ^ ' 

The exact expressions for the extinction probabilities are 

A = ^, 



pjw = <! i^^^nw. (1-6) 

A-Aic-HA-P)' 

From (1.2), we can easily infer that the probability generating function of Ny{t\t > 0, 

G^(u,t) =Ew^-(*', |it|<l,0<i^<l,i>0, (1.7) 
satisfies the Cauchy problem 

d" d 

— G^{u,t) = {\u- ^i){u~l)—GAu,t), i^e (0,l],|u| < 1, 
G^{u, 0) = u. 

We will show below that from (1.8), one can arrive at the subordination relationship 

N,it)'^- N{T2,it)), t>0, (1.9) 

where T2^{t), t>0, is the random time process whose distribution is obtained by folding 
the solution of the following fractional diffusion equation: 
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^qix,0)=Six). 



0<i^< 1, x-eR, i>0, 



(1.10) 



The process N{t),t > 0, found in (1.9), is the classical hnear birth-death process whose 
distribution is given in (1.4), (1.5) and (1.6). A relationship similar to (1.9) also holds for 
the fractional pure birth process [8] and the fractional Poisson process [2]. In this context, 
it represents the main tool of our analysis and leads to a number of interesting explicit 
distributions. We consider the subordinator related to (1.10) because the probability 
generating function of the distribution of (1.9) satisfies the simplest fractional equation 
generalizing the classical one. 

For the extinction probabilities of the fractional linear birth-death process, we have 
the following attractive formulas: 



/i A — /i 



+ 00 

E 



Poit)={ 



1 - 



/i — A 



^.,i(-i''(M-A)m), A<M, 



(1.11) 



-f OO 



c-'"E^,i{~Xfw)dw, 



for t >0,0 < 1. 

The function Ea^p{x) appearing in (1.11) is the generalized Mittag-Leffler function, 
defined as 



+ 00 



r(am + /3)' 



xgR,q:>0,/3>0. 



(1.12) 



From (1.11), we can easily retrieve the classical extinction probabilities (1.6) for v ~1 
by keeping in mind that £"1^1(0;) — e^. 

For the state distributions p^{t)^ t > O, k> 1, we have formulas similar to (1.11), but 
with a more complicated structure: 



A — /i 



E 



l + k 
I 



I fe-i 



E(-ir 



r=0 



k-1 

r 



Pkit) = <^ 



fc-i / > \ 2 00 
/i — A 



E 



l + k 
I 

k-l 



xE,,i{-{l + r+l){X-fi)e), 
k-l' 



A>/i, 



(1.13) 



xE,,ii-il + r + l){^i~X)t''), 



X<n, 
X — 11. 
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Also from (1.13), for v = 1, one can reobtain the distributions (1.4) and (1.5). 

We will show below that the probabilities p'^{t), t>Q, > 1, appearing in (1.13) are 
strictly related to the distributions of the fractional linear pure birth process M^{t), t>0, 
with an arbitrary number of progenitors and a birth rate equal to A — /i with A > /i. In 
particular, we can extract from the first line of (1-13) that 



Pi{N,{t) = k\N,{0) = l} 



OO r 



A ^ 

(=0 



(1.14) 



X Pr{M^(t) = k + p/„(0) = / + 1}, 



where 



is a geometric law for the number of progenitors. We also note that for A = /i, the 
distribution (1.13) can be expressed in terms of the extinction probability (1-11) by 
means of 

FT{N,it) = k} = ^^^^^^^^-^[X{l-Pr{N,{t)=0})l fc>l,i>0. (1.16) 

The extinction probability (1-11) can be viewed as being a suitable weighted mean of the 
waiting times of the fractional Poisson process P^{t),t > 0, for which it is well known 
that [2] 

Pr{A/'7>.^(i)=0} = £;,,i(-Ar), t>0. (1.17) 

The fractional linear birth-death process dealt with in this paper provides a general- 
ization of the classical linear birth-death process and may well prove to be capable of 
modeling queues in service systems, epidemics and the evolution of populations under 
accelerating conditions. The introduction of the fractional derivative furnishes the sys- 
tem with a global memory. Furthermore, the qualitative features illustrated in the last 
section show that the fractional counterpart of the linear birth-death process has a faster 
mean evolution (and variance expansion), as was pointed out in similar fractional gener- 
alizations, for example, for the Poisson process (see [2, 3, 5, 11]), for fractional branching 
processes [10] and for pure birth processes [8]. 



2. The extinction probability of the fractional linear 
birth— death process 

We begin this section by proving the subordination relationship (1.9) which is relevant 
to all of the distributional results of this paper. 
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Theorem 2.1. The fractional linear birth-death process N,y{t),t > 0, can be represented 



N,{t)^N{T2u{t)), t>Q,0<v<l, 



(2.1) 



where N(t), t> 0, is the classical linear birth-death process and T2v{t),t > 0, is a ran- 
dom process whose one- dimensional distribution coincides with the folded solution of the 
fractional diffusion equation 



q2v q2 



(2.2) 



subject to the initial conditions g(x,0) = 5{x) for Q <1 and also qt{x,0) = for 1/2 < 
v<l. 

Proof. The Laplace transform Gi,{u,z) = e~^*G^{u,t) dt, applied to the fractional 
PDE 



yields 



|-G,(u,t) = [Xu - ii)iu ~ l)-^G4u,t), < < 1, 
Gv{u, 0) = u, 



d 



(2.3) 



z''G^{u,z)-z''-^u^{Xu~fi){u^l) — G^{u,z), < < 1, z > 0, |u| < 1. (2.4) 



du 



We now observe that 



G^{u,z) : 



'^u''PT{N^{t):=k} 



Lk=0 



dt. 



(2.5) 



If (2.1) holds, then 
G^(u,z) 



e-^'* / Pr{Af(s) = k) Vv{T2u{t) € ds) 

.k=o -^0 

) nOO 

e-^* / G{u,s)Pr{T2,{t)eds} 
Jo 

G{u,s)z''-h-'''' ds. 



dt 



(2.6) 



In the last step, we applied the folded version of equation (3.3) in [6] for c—1, that being 
therefore 



f 

Jo 



e-^*Pr{T2^(i)eds} = e"''" z^^'ds 



(2.7) 
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We now show that (2.6) satisfies equation (2.4); by inserting the Laplace transform into 
(2.4), we obtain 



G(u,s)e-''- ds- z''-^u = {\u- ii){u-l)z''-^ / _G'(w,s)e~"" ds. (2.. 

Jo ou 

The inversion of the integral with d/du is justified because 

oo 

Y2ku^'^Vi{N{s) = k} 



k=l 



(2.9) 



< ^ fc Pr{iV(s) = k} EiV(s) < cx). 



fe=i 



Taking into account that G{u, t) satisfies the first-order PDE 



dG dG 

97 ^(^"-'^^("-^^a^' 



(2.10) 



from (2.8), wc have that 



G(w,s)c '^^ ds-w= / —G(u,s)e ''^ ds 
Jo OS 



Giu,s)c-''"\lzT+z" I G{u,s)c-''^ ds (2.11) 



G'(M,s)e-'*^"'ds. 



This shows that (2.1) holds for the one-dimensional distributions. This concludes the 
proof of Theorem 2.1. □ 

Remark 2.1. For = l/2",n e N, the density /^a,^ of the random time T21, appearing 
in (2.1) becomes the probability density of an [n — l)-itcratcd Brownian motion, that is, 

Pr{ri/(2„-i)(i) e ds} = Pr{|Bi(|S2(. • • \B,M ■ ■ -m e ds} 



°° p-s /(4wi) 



^/AlXU^l 

\/47rt 



dioi 



■du2--- (2.12) 



dw„-i. 



as can easily be inferred from [7], Theorem 2.1. The difference between (2.12) and its 
corresponding formula in the cited paper is that here, the diffusion coefficient is equal to 
one. 
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In the following theorems, we separately derive the three different expressions of the 
probability of extinction in the cases X > fi, X < fi and A = /i. Wc prefer to treat them 
separately because their proofs are somewhat different. 



Theorem 2.2. For a fractional linear birth-death process with rates X> fi, the proba- 
bility of extinction has the form 



p[;(t)=Pr{iV,(i)=0} 



E 



(2.13) 



A A 



for t > 0,0 < V < 1, and where E^_i{x) is the Mittag-Lejfler function (1.12). 



Proof. In light of the subordination relationship (2.1) of Theorem 2.1, and by taking 
into account the extinction probability of the classical linear birth-death process 



Pr{iV(i) = 0} = 



/i — /ze 



-t(\-^l) 



A - /xe-*(^-'^) ' 



wc can write that 



Pr{iV,(t) = 0} = 



t>0, 



Vr{T2u{t)(^As} 



(2.14) 



for alH > and < < 1. By taking the Laplace transform of (2.15), we obtain that 



(2.15) 



/ e-^*Pr{7V^(<) = 0}dt 
Jo 



X - /ie-'*(^-'") ' 



A7o 



oo 

(l-e-"(^-^)) 

rn—0 

oo 



ds 



(2.16) 



m=0 



.,1^-1 



^Q-s(\-n)m-sz" _ Q-s{\-^l){l+m)-sz" 
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The above steps are valid because < ^e ^'^^ < 1 for A > /i. By keeping in mind the 
Laplace transform of the Mittag-LefHer function E^^i{—xt'^), 



e-''*E^^i{-xt'')dt: 



we readily arrive at the claimed result. 



(2.17) 
□ 



Remark 2.2. When i/ = 1, we obtain from (2.13) the form of the extinction probability 
(2.14) for the classical birth-death model: 



Pr{iV(t) = 0} = 



/.t — A 
A 

fi — X 
A 

/i — A 
A 

/i — ^ic 



-foe 

E 

7n—l 



1 



^— (A — 



- 1 



A 



1 - (^/A)c-(^-")* 

(A^/A)e-*(^-^) 
1 - (Ai/A)e-*(^-A') 

-t(A-A() 



A 



(2.18) 



A 



A - ^c-*(-^-^) ■ 

From (2.15) for v =1, Pr{T'2(t) G ds} = 5{s — t) and we again retrieve result (2.14). 
Remark 2.3. From (2.13), we note that 



PT{N,{t)=0} 



^ A 



Vi.e(0,l], 



(2.19) 



which is the asymptotic extinction probability, irrespective of the value of v. 



Let us now deal with the case A < /x, that is, when the rate of birth is strictly lower 
than the rate of death. 

Theorem 2.3. For /i > A, the probability p^it) — PT{N^{t) — 0} of complete extinction 
of the population is 



/i — A 



m=l 



(2.20) 



where t> 0, 0<i'<l o,nd E,j^i{x) is the Mittag-Leffler function (1.12). 
Proof. We start by rewriting (2.14) as 



(2.21) 
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Using (2.1), we are able to write 



Pr{r2,(t) e ds}. 



By applying the Laplace transform to (2.22), wc obtain that 



z^-h-"""" ds 



(A/^)e-«(A'-^) - 1 

•^0 rn=0 



-s(/j-A) 



+ 00 



m=0 

+ 00 



m /• + CX3 



(1 - c-^(A'-A))g-s^''g-s(M-A)m 



m=0 

+ CXJ 



~s(/i— A)7Ti— sz" _ g-s(/^— A)(m+1) — S2"^^ 



1 



1 



l,m=0 

I, m=l 



i yfi 

.Vti -I- rl' .^-^ \ // 



1 



+ 00 



1 



\ u I [u — X)m + z'^ A •■'-^ V /i / (u — X)m + 



E 



1 



(2.22) 



(2.23) 



1 - 



+ 00 

E 



E 



1 



r'^VM/ (a« - A)m + 



Inverting (2.23) by means of (2.17), we retrieve formula (2.20). 



□ 



Remark 2.4- When v = 1, we reobtain from (2.20) the extinction probability of the 
classical birth-death process (2.21): 



Plit) = 1 



+ 00 



- {i.i — X)mt 
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1 
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(2.24) 



/i — A 



(A/a^)c- 



-\{p-X)t 



1 - (A/^)e-(^-^)* 



= 1 



(A2^)e~^^"^^* - Ae-(^-^)* _ A - Ae-(f-^)* 



A - (A2//^)e-(^-^)* 



A-(A2/M)e-(^-^)* 



I _ e-(A'-->')t 



1 - (A/Ai)c-(A'-^)* Ac-*(^'-^) - n ■ 

Remark 2.5. Population extinction in the long run is evident from (2.20) as 

Pr{iV,(i)=0}*^°°l, (2.25) 
due to the death rate exceeding the birth rate for all < < 1 . 

In the next theorem, we treat the remaining case, that is, when /i = A. 

Theorem 2.4. For the fractional linear birth process, when the rate of birth equals the 
rate of death (i.e., when X — fi), the extinction probability pQ{t) reads 



\-iv r+oo n + oo 

Po(i) = / e-'"E^,^{-wXf)dw^l- c-'''E^,i{-Xt"w)dw 

^ Jo "'0 

with t> 0, 0<i'<l and where Ei,^i{x) is the Mittag-Leffler function (1.12). 
Proof. Again using (2.1), we write 



(2.26) 



Po(i) 



As 



■Pr{T2,(t)eds}. 



1 + As 

We now apply the Laplace transform once again, thus obtaining 

le- 



(2.27) 



As + 1 

se^'' 



■ds 







u-l 



Xz 



(t/=s(z^+Aiu)) 

— AZ 



+ 00 



+ 00 



+ 00 



+ CJO 



- s — wXs 



ds dw 



(2.28) 



y 



+ 00 



1 



z"" + Xw z"" + Xw 



z^ + Aw z'^ + Xw 



■ dw. 



■ dw 
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By inverting the Laplace transform, we obtain the integral form 

^+00 



p^{t) = \ e-"" u''-^E^^^{^wXu'')E^^i{^wX{t-uY)dudw, (2.29) 
Jo Jo ' 

which involves convolutions of generalized Mittag-Leffler functions defined, for 

example, in [9], equation (1.56), page 17. The inner integral in (2.29) can be worked out 
explicitly as follows: 

t 

u''-^E^^^{~w\u'')E^^i{-w\{t - uy)du 

^ ^T{vm + v)T{vr + l) Jo ^ " 

m=Or=0 ^ ' J \ ' ) JO 

^ - - {-n,X)m+r ^^^^^^^^^ T{vm + v)T{vr + l) 

^ ^T(um + v)T(vr + l) T(iy(m + r) + ly + 1) ^' ' 

m— r— / \ / \ \ / / 

0000 y -.xyj 00 / .-.^ 

^^T(un^u^l) ^r(i/n + z/ + l) ^ ^ 

m—On—m ^ ' n— ' 

The extinction probability now reads 

p'^it)^ / e-'"E^^^{-wXf')dw. (2.31) 

^ Jo 

Using the relationship 

^E,^i{x)^-E,,,{x), (2.32) 
da; v 

the extinction probability (2.31) takes the alternative form (2.26) because 

Poit) -— rE^yr^^'^-^dy 

^ Jo 

= - f i?...(2/)e^/(^*")d2/'^^^^ r c«/(^*'^)A£;,,i(y)dy (2.33) 
= r e^/(^*")i?,i(y)dy^"'=-?L/(^*''»l- ro-E^A-^t'^w)dw. 

J-oo Jo 

This completes the proof of (2.26). □ 
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Remark 2.6. From (2.26), when v = 1, we again retrieve the classical form 

.oW = ^, (2.34) 

as expected. 

Remark 2. 7. The limiting extinction probability when = A is 

Pr{iV,(t) = 0}*^°°l (2.35) 

for all values of < < 1 . 



Remark 2.8. The last expression in (2.26) is in some ways similar to the Riemann limit 
for ^ ^ A of (2.13) and (2.20). 

Remark 2.9. We can rewrite the probabilities (1.11) in an alternative form which per- 
mits us to give an interesting interpretation to their structure. 
For the case A > /U, we can write 



^ A 



1 - J2 MS = m\g > IjE^^ii-t^'iX - n)m) 



(2.36) 



where Q is a geometric r.v. with distribution 

Pr{Q ~ m) A — /i/ /i\™A 



Pr(g = m\g > 1) 



Pr(g>l) 



A V A 



m > 1. 



(2.37) 



The treatment of the opposite case X < fi is similar except that a different conditional 
geometric r.v. Q' must be introduced, defined as 



Pr(^'^H^'>l) = ^(^)'"^, 



m> 1, 



and thus 



(2.38) 



(2.39) 



A well-known property for a fractional Poisson process J\f^{t),t > 0, of degree < < 1 
and parameter A > is that [2] 



Pr{AC(i) - 0} = £;,,i(-r A) = Pr(r. > t), 



(2.40) 
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where Tu = inf(s : A/'i/(s) = 1)- This permits us to rewrite the extinction probabihties also 
in terms of waiting times of a fractional Poisson process with a random rate XQ . 

For the case A = /i, the interpretation is straightforward because the waiting time of the 
related fractional Poisson process has a rate X£, where £ is an exponentially distributed 
r.v. with parameter equal to one. 



Remark 2.10. In the case ^ = A, it is well known that the extinction probability in 
the classical birth-death process, Pq{s)^ s > 0, satisfies the nonlinear Riccati differential 
equation 

p'a[s) + 2Xpo{s) = A + X[po[s)f. (2.41) 

By using (2.41), we can provide an alternative proof for the subordination relationship 
(2.1): 



p;,(s)Pr(r2.(t)eds)+2Ap^(t) = A + A / [po{s)YVr{T2.{t) e ds) 



A 



(1 + A.s)2 

■ (1-AV) 
(1 + As)2 

1 - As 



Pr(T2^(t) e As) + 2Xpl{t) = A + A 
Pr(T2,(t)ed.s) = A-2Ap^(t) 



(1 + As)- 



(2.42) 



■Pr(r2,(t)eds) 



1 + A.S 
2Po(i)-2 



Pr(r2,(t)eds) = l-2p|;(i) 
°° Xs 



1 + As 
As 



Pr(r2.(t)eds) 



1 + As 



Pr(r2,(t)Gds). 



(2.43) 



Remark 2.11. By exploiting the subordination relationship (2.1) and the fact that the 
extinction probability in the classical case satisfies the integral equation 



Po{t) 



(2.44) 



we can give an integral form for Po(i): 



PoW = 



+ 00 



(^+^)"/idu+ / Ae-(^+^)"[po(s-it)]^dulPr{T2^(i)eds}. (2.45) 



(2.46) 



We note that the first integral of (2.45) can be worked out explicitly as follows: 



"'0 



Pr{T2^(t) edsjdi 



1 



z A + /i + 
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This can be directly inverted so as to obtain 

1*00 j*S 

/ / e-(^+^)"/iPr{r2^(t) edsjdM 
Jo Jo 

(2.47) 

(_(A + Ai)t-)™ 



m=0 



A + /i 



(to + i)r(j/TO + v) 
[i-i;,,i(-(A + /i)t'^)]. 



3. The state probabilities of the fractional linear 
birth— death process 

Here, we present three theorems concerning the structure of the state probabihties 
Pr{A^i^(t) = fc}, i > 0, with < < 1. Three cases must be distinguished and treated 
separately, as in Section 2, namely A > /i, A < /i and A = /i. 

Theorem 3.1. For the case A > /i, the state probabilities p'^{t), k>l, t>0, 0<!^<1, 
in the fractional linear birth-death process N^{t), t>Q, have the following form: 

(3.1) 

X £;,,!(-(/ + r + l)(A- 

Proof. By exploiting the subordination relationship (2.1) and conveniently rewriting the 
well-known form of the state probabilities of the classical linear birth-death process, we 
have that 

Pm^iX-i.fx'^-' --''-'' \\_^,-,.Js),^^ Pr(r,.(t)ed.). (3.2) 

By applying the Laplace transform, we obtain 

poo fl _ p-(A-A')s\fc-l 

^ y (A- Aie-(-^-^)^)'=+i 



7— n r—n V / \ / \ 



i=0 r=0 



r 



/•CXD 

/ e-''^"'e-(^-'')^'e-(^-'')"''e-(^-^'"ds (3.3) 
Jo 
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which can be easily inverted by using (2.17), thus obtaining (3.1). □ 

Remark 3.1. We check that for v formula (3.1) converts into the well-known dis- 
tribution of the linear birth-death process, thus being its fractional extension. For 1^ = 1, 
we get from (3.1) that 

pi^^={^)'f:{'^j''){^Tf.^-^r{^-^^ (3.4) 



A 

We now observe that 

fc-i 



/=o ^ / \ / ,.^0 



^(-1)'' ( / ) e-^f^-^)-^ = (1 - e-(^-^)*)''"\ (3.5) 

oo / , , \ / \ i / \ — (fc+1) 

where, in (3.6), we applied the binomial expression 

OO / \ GO / \ 



This permits us to write 



A ; (l-(/i/A)c-(^-f)t)'^+ 
which coincides with (1.4). 

Remark 3.2. In order to prove that Yl^=iiPki^) ~ 1 ^O'' A > /i (formula (3.1)), we can 
again apply the Laplace transform and prove that X^fe^o /o°° *^^^*^'fc(^) = l/z. We first 
calculate 



, \ 2 OO fe— 1 , 7 \ / , -, \ / \ ' 

A — M \ I + k\ k — 1\ , -.xr/M 



fc=l ^ ^ l=Q r=0 
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fc=l \ / ;=o r=0 



By keeping in mind formulas (3.5) and (3.6), we have that 

oo poo 

/ o~^Vk{t)At 



fc=i 
(A-^)z''-i 



A - /ie-'^t^-'^) ■ 



By using the Laplace transform of the extinction probability (second line of formula 
(2.16)), we finally obtain 



c-^Yk{t)dt^{X-f,)z'^-^ 



A - ^je-'*(^-A') 







A - ^e--'*(^-A') 
1 

z^'-ie""^ ds=-, 
z 



.s(A-p) 

z'^-^e""^" ds (3.11) 



as desired. 

Remark 3.3. The distribution (3.1) can be expressed in terms of the probability law of 
a fractional linear birth process with rate X — (J., which reads 

q'^it) = PT{M,At) = fc + l\AUO) =1 + 1} 

where Z + 1 initial progenitors are assumed (see [8], formula (3.59)). If we write 

\ / \ ' 



Pr{G = ;}=(l-Jl( J) , />0, (3.13) 
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then formula (3.1) can be rewritten as 

Pkit) 



A — /i ^ ^-v I + k f /i 



A 



A 



E 

(=0 



k V A 



PT{AUit) = fc + /|A/^.(0) = / + 1} 



Pr(G = + T3^Pr(G = /) 



k{X-ii) 

X Pr{A/i.(i) = fc + p/^(0) = / + 1} 



(3.14) 



because 



\J\\J k{\-^i) 



Pr(G = /). 



(3.15) 



Result (3.14) shows that for large values of k, we have the following, interesting, approx- 
imation: 



X — fi 



J2 Pr(G = Pr{A/.(i) - fc + l|A/,(0) = 1 + 1}. 



1=0 



Theorem 3.2. For a fractional linear birth-death process N^{t), t>0, fJ.> X, the prob- 
abilities p'j^(t) = PT{N^{t) = k}, k >1, have the following form: 



Pkit) = 



/i — A \ fx 



fc— 1 oo 



E 

1=0 



l+k\ fx 



fc-1 



E(-i) 



r=0 



r fk-l 

r 



(3.16) 

E,,,{-{l + r + l){fi-Xy). 



Proof. By again using relationship (2.1), thanks to formula (1.4) suitably rearranged, 
we can write 

roo /'„-(/j-A)s _ 1 ^|fe-l 

P^t) =1 - ^)'^'^'"^^^'" (tal^3073^ W e d.s). (3.17) 

By applying the Laplace transform, we have (omitting here some steps similar to those 
of the proof of the previous theorem) 



(-^)'=+i(l- (A/Ai)e-(A'-^)-^) 



k+l 



j^-ie-"^" ds 



A\VA 



fc— 1 CO 



E 

i=0 



i+k\ fx 



I fe-1 



E(-i)'- 



, /fc- 1 



r=0 



(3.18) 
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By transforming equation (3.18), we easily arrive at the result (3.16). □ 
Remark 3.4- When k ~ 1, equation (3.1) takes a simple form: 



■. \ 2 oo / \ 1—1 



(3.19) 



1=1 



where A > /i. For the ease A < /i, we obtain essentially the same expression with A and fi 
exchanged. 

An interpretation similar to that in (3.14) is valid for the case ^ > A as well. The 
following theorem describes the structure of the state probabilities p\{t), > 1 , in the 
case where ^ X, that is, when the birth rate equals the death rate. 

Theorem 3.3. In the case /i = A, the probabilities p'^{t) ~ Y'T{N^{t) = fc} of the frac- 
tional linear birth-death process read 

Vv{NS) = fc} = ^^^^1T^^[^(1 (3.20) 
with fc > 1 and i > 0. 

Proof. The explicit form of the distribution Pr{iV^(t) = > 1, for the fractional 
linear birth-death process, in the case A = /z, can be evaluated in the following manner. 
In light of (1.9), we have 

/>oo 

Pr{7V,(f) = fc}= / Pr{iV(s) = fc}Pr{T2,(i)eds} (3.21) 
Jo 



so that 



poo 

Ll{z)^ / e-^*Pr{7V^(i) = /s}dt 

"'0 



(As)'' ^ sz" 



(3.22) 



(1 + As) 



fc+i 



z^-'e-"'' ds. 
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This is because for the \ = n case of the classical birth-death process, we have that (see 
[1], formula (8.53), page 95) 



Pr{iV(t) = k] = 



iXt) 



fc-i 



(l + At)^+i' 



fc> 1. 



(3.23) 



We note that the extinction probability cannot be extracted from the above formula since 
it reads 



Fi-{N{t) = 0} 



Xt 



l + Xt 



(3.24) 



This implies that we have a different expression for A: > 1 and fc = for the fractional 
linear birth-death process as well. 

Formula (3.22) can be expanded out as 



{- 




-1 d*^- 


nOO 






kl 


dXf^J 







i- 




-1 ^fe 








kl 


dAfc j 


[ 




(- 


-1)'=A'=- 


-1 ^fc r 




1 




k\ 


dA'^- 


Jo 


Jo 


(- 












kl 


dA^ 





1 



ds 



1 A 



s 1 + As 



(e'^'"'- - Ae-'"(i+^-^))z'^-ie~""'' dsdw 



(3.25) 



/q wX + , 

By inverting the Laplace transform, we have that 



■ dw 



iw 



w + z" 



Pi{N,{t)^k} = 



{- 




^k 






k\ 


dA*^- 


IL 


i- 


_l)k-lxk 


-1 ^k 




kl 


dA'^ 


{' 


,l)k-lxk 


-1 d^ 




kl 


dA^^ 



{E^^ii-wt'') ~ Xe-"" E^^i{~Xwt''))dw 
X / e-"'£;^,i(-Awi^)dw 



Formula (3.20) is thus proved. 



(3.26) 



□ 



It is important to note how all the state probabilities p^(i) depend on the extinction 
probability Po(t). 

Remark 3.5. For = 1, we can extract from (3.20) the classical formula (3.23) because 



pl{t)^Vr{N{t)^k}^ 



k—1 \ fc— 1 jfc 



(-l)'=-iA 



k\ dA'^ 



A 

l + Xt 



(3.27) 
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and because 



dfc r 



A 



l + \t 



E 

J=0 



A 



j / dA^ dA*:--? \ l + \t 



1 



A- 



d* 



1 



jfe-i 



lA''^ V 1+At / dA*^-! \l + \t 



1 



, (-lfk\t^ , {-lf-^{k-l)\t^-^ 
A— + k- 



(l + At)'=+i 



(1 + At)fc 



f;+'i'P (-l)''M-Afct+(l + At)fc] 



(1 + At)'^+ 



-(-1) 



fc-i 



(3.28) 



Remark 3.6. From the representation on the last hne of (3.26), it is possible to give an 
alternative proof of the subordination relationship (2.1) when A: > 1, as follows: 



k—1 1 Ak 



kl dA*^- 
fc! dA^ 



X- 



(3.29) 



fc! 



Exploiting (3.28), we readily obtain 



°° d^ 
dA^ 



1 + As 
A 



Pr(r2.(t)eds) 



1 + As 



Pr(r2,(0eds) 



(3.30) 



Pl(t) 



(As)^- 



(1 + As)fe+ 



rPr(T2.(i)ed.s). 



(3.31) 



Remark 3. 7. Here, we present two other interesting relationships. The first one is simply 
a particular case of formula (3.20) when fc = 1, that is, the probability of having one 
individual in the process at time t is 



PT{N,{t) = i} = ^[x{i-p^,m 



(3.32) 



The second relationship is again a particular case of formula (3.20) with u = 1/2. In that 
case, recalling that 



(3.33) 
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we obtain 



Pr{iVi/2(t) = fc} = ^^^^ _[A(1 

_ (_l)fc-iAfe-i 
^ fc! 



V^Ja Jo 



fc! dA*' 
fc! dA*^ 
fc! dA*^ 



2A 



^/7T 7o 1 + 2Ay^/^ 



dy 



2A / ^ ^d«; 



1 + \V2w \/27rt 
2A 



1 + XV2B{t) 



(3.34) 



where B{t), ^ > 0, is a standard Brownian motion. 



4. Some further properties of the fractional Unear 
birth— death process 

The analysis of the moments of the fractional linear birth-death process gives us useful 
information concerning the behaviour of the system. Starting from (1.8), we easily see 
that 



(4,1) 

11=1 



is the solution to 




By again applying the Laplace transform, we have that the solution to (4.2) is 

EN,{t)^E,s{{X-fi)n, t>0. (4.3) 

In the case X> fi, the result (4.3) shows that the mean size of the population coincides 
with that of a fractional linear pure birth process with rate A — /i > (see [8]). Result 
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fe=0 



¥.N^{t) = ^kVT{N^{t) = k} 

fe=0 

^ fc / Pr{iV(s) = k} Pr{r2,(<) e ds} 

7,,_n J 

3(A-M)3pr{r2,(t)eds}. 

The Laplace transform of (4.4) yields 

/•OO 

*^N^{t)At^ / e(^-^)"z'^-ie-""''ds 



(4.4) 



(4.5) 



and this confirms (4.3). 

By again applying (1.8), it is also possible to derive the variance Wa.rNi,{t), t > 0, 
of the number of individuals in the population at time t. We start by evaluating the 
second-order factorial moment /i(2)(i) = ^[Ni,{t){Ni,{t) — 1)], t > 0. From (1.8), after 
some straightforward steps, we see that 



^(2)(t)=E[iV,(t)(7V,(i)-l)] 



is the solution to the following differential equation: 



— Ai(2)(i) = 2AE7V,(t) + 2(A - Ai)/^(2)(0, < < 1, 

,/^(2)(0)-0. 



(4.6) 



(4.7) 



In order to solve (4.7), we apply the Laplace transform, obtaining, in the case A 7^ /x, 

..^-1 1 



>(2)(t) At = 2A 7- — 

' z" -{X- ^i) - 2(A - ^) 



2Az 



A 



z" - 2(A - /i) z'' - (A - /i) 



(4.8) 



The Laplace transform (4.8) can be inverted, thus leading to the explicit expression of 
the second-order factorial moment as 



M(2)(0 



2\ 
A — ^ 



[S,,i(2(A-M)t")-^..i((A-/^)i")]. 



(4.9) 
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From the first expression of tfie Laplace transform in (4.8), we also have that 

M(2)(t)=2A / s''-'E,^,{2{X- ^i)s•')E,,li{\- ^i)it- sy)ds. (4.10) 





By applying similar calculations to those of (2.30), we prove result (4.9). 
From (4.9), we can easily write that 

2A 

VariV,(t) = - [E,,,i2{X- - E,^i{{X- fi)f')] 

A — /i 

+i?.4((A-/.)n-<i((A-/^)n (4.11) 



i?.,i(2(A - ^l)n - ^i?.,i((A - ^^)n - Et,,{{\ - 



X — ^ ' X — II 



Remark 4-1- When ;y = 1, wc obtain from (4.11) the expression for the variance of the 
classical linear birth-death process as follows: 

VariV(t) - ^e2*(^-^) - A±^e*(^-^) - e''^^->'^ 
A — u X — 11 

. (4.12) 

X — f-i X — /-/ 

Remark 4-2. When = 0, that is, in the case of pure linear birth, we obtain from (4.12) 
the expression of the variance of the classical linear pure birth process and from (4.11) 
that of the fractional linear birth process (see [8]). 

In the case A = /i, from (4.7), it is easy to show that 

2Xf'^ 

Var7V,(t) = j,^^, (4.13) 

in accordance with the well-known result of the classical linear birth-death process for 
X = n which reads Var7V(t) = 2At. 

Remark 4-3. Wc can directly evaluate the mean value EN{t) ior X = ^ in the following 
way: 



(4J4) 

z=Xt/{l + Xt) 



-(')^g<(I^)^^g<T, 



(1 + At)2dz 

k—1 



(1 + Ai)2 dz 1 - z 

z=At/(l + At) 
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1 1 

" (1+At)2 (l-z)2 

The assumption that \ = ji imphcs that the mean size of the population 'EN^^t), i > 0, 
is equal to one (number of original progenitors) for alH > and all < < 1 (this is also 
confirmed for A = /i by (4.3)). 
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